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The dynamical spin response of doped two-leg ladder antiferromagnets is investigated based on 

■ the fermion-spin approach. Our calculations clearly demonstrate a crossover from the incommensu- 
rate antiferromagnetism in the weak interchain coupling regime to commensurate spin fluctuation 

^ in the strong interchain coupling regime. In particular, the nuclear spin-lattice relaxation rate 

extracted from the commensurate spin fluctuation decreases exponentially with decreasing temper- 
atures. The behaviors of the spin dynamics in the strong coupling regime are quantitatively close 
. to the experimental results of Sri4_ :E Ca :E Cu2404i . 
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In recent years the novel two-leg ladder antiferromagnet Sri4Cu2404i, being situated between doped one- and two- 
dimensional antiferromagnets, have been experimentally investigated as well as theoretically 1 . This followed from the 
fact that this material has a spin liquid ground state 2 , which may play a crucial role in the superconductivity of doped 
cuprates as emphasized by Anderson 3 . Once carriers are added to the material Sri4Cu2404i, such as the isovalent 
substitution of Ca for Sr, a metal-insulator transition has been observed 1, 4,s , and further, this doped two-leg ladder 
material Sri4_ x Ca2;Cu2404i is a superconductor under pressure in low temperatures 4,5 . All cuprate superconductors 
O \ found up now contain square Cu02 planes, whereas Sri4_ :E Ca ;E Cu2404i consists of two-leg ladders of other Cu ions 
and edge-sharing Cu02 chains 1,4 ' 5 , and is the only known superconducting copper oxide without a square lattice. 
| Experimentally, it has been shown by virtue of transport measurements that there are the regions of parameter space 
. where the resistivity is linear with temperatures 5 , one of the hallmarks of the exotic normal state properties found in the 
two-dimensional cuprates 6 . However, NMR and nuclear quadrupole resonance (NQR), particularly inelastic neutron 
scattering measurements in the same regions of parameter space indicate that Sri4_. E Ca a ;Cu2404i in the normal state 
is an antiferromagnet with the commensurate short-range order 1,7,8 . This commensurate spin correlation is energy 
dependent, and persists in a wide range of doping 1,7,8 . Moreover, NMR and NQR spin-lattice relaxation rate extracted 
■ from this commensurate spin fluctuation decreases exponentially with decreasing temperatures 7,8 . These magnetic 
behaviors are different from these of the doped two-dimensional high-T c cuprates 9,10 , where the incommensurate 
spin fluctuation for the single layer lanthanum cuprate 9 , and both low-energy incommensurate spin fluctuation and 
high-energy commensurate [tt,tt] resonance for the bilayer yttrium cuprate 10 in the normal state are observed. 

Theoretically there is a general consensus that the unusual physical properties of the two-leg ladder materials 
are due to the quantum interference between the chains in the ladders 1 . Applying the bosonization procedure to 
the two-leg t-J and Hubbard ladders, it is shown 11 that all spin excitations are gapful and a singlet pairing be- 
comes the dominant instability. These results are confirmed by some research groups within the different theoretical 
frameworks 12 . Furthermore, based on the antiferromagnetic Heisenberg model, the susceptibility and spin-lattice re- 
laxation rate have been discussed 13 , and the results show that the large contribution to the spin-lattice relaxation rate 
comes from processes with the wave vectors around the antiferromagnetic zone center. Within the one dimensional 
gapped quantum nonlinear a model, an effective classical model has been developed to study the spin transport of 
the two-leg ladder antiferromagnets 14 , and the result obtained for the spin-lattice relaxation rate is close to what is 
experimentally measured. However, to the best of our knowledge, no systematic calculations have been performed to 
show why the commensurate neutron scattering peak can be observed in the doped two-leg ladder antiferromagnet 
Sri4_ x Ca^Cu2404i? This is a challenge issue since it is closely related to the doped Mott insulating state that forms 
the basis for the superconductivity 4 . In this paper, using the fermion-spin theory 15 which implements properly the 
local single occupancy constraint, we calculate explicitly the dynamical spin structure factor within a, t-J ladder and 
show that in the regions of parameter space given by the experiments, one can reproduce all main magnetic features 
observed experimentally on Sri4_ ;E Ca :c Cu2404i in the normal state, including the energy dependence of the neutron 
scattering peak position and exponential decrease of the nuclear spin-lattice relaxation rate. 

The two- leg ladder is defined as two parallel chains of ions, with bonds among them such that the interchain 
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coupling is comparable in strength to the couplings along the chains, while the coupling between the two chains that 
participates in this structure is through rungs 1,4 . The essential properties of the doped two-leg ladder antiferromagnet 
can be described by the t-J ladder as 16 , 

H = -h E cLc i+f)aa t ± E(cL^2 CT + h. c .) 

ifjacr icr 

- ^E CiggCiaa + J|| ^ S ia ■ S i+ fj a + J± ^ Sjl • S i2 , (1) 

iaa vqa i 

with the local constraint C\ aa Ci aa < 1 to remove double occupancy of any site, where r) = ±cox, Co is the lattice 
constant of the two-leg ladder lattice, which is set as the unit hereafter, i runs over all rungs, <r(=T, I) and a(= 1, 2) 
are spin and leg indices, respectively, C\ aa (Ci arT ) are the electron creation (annihilation) operators, Sj = C\ a uCi a j2 
are the spin operators with a = (p x ,Oy,o z ) as the Pauli matrices, and /i is the chemical potential. For the two-leg 
ladder materials, it has been shown from experiments 1,7 ' 8 that the exchange coupling J|| along the legs is greater 
than exchange coupling J±_ across a rung, i.e., Jm > Jj_, and similarly the hopping t|| along the legs is greater than 
the rung hopping strength t±, i.e., tu > tj_. On the other hand, the strong electron correlation in the t-J model 
is reflected by the local constraint 3 , which can be treated properly within the fcrmion-spin theory 15 , where electron 
operators Ci a \ — h\ a S^ a and Ci a [ = h\ a S^ a are represented by the spinless fermion operator hi a carrying the charge 
(holon) and the pscudospin operator Si a representing the spin (spinon) , then it naturally incorporates the physics of 
charge-spin separation. In this case, the low-energy behavior of the t-J ladder (1) can be rewritten in the fermion-spin 
representation as, 

H = H E h Uahia{St a Sr + ~ a + Sr a S+ +rla ) + 1 ± E(4fc2 + h\ 2 ha)(S+S. 2 + 
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h\ a h ia + J|| off &ia • S i+ jj a + J± cS Sil ' Sj2, (2) 



where J|| e ff = J\\ [(1— p) 2 — </>|], Jj_ e fi = </[(l— p) 2 — </>j_], P is the hole doping concentration, the holon particle-hole order 

parameters <f>u = (/ij a /ij+^ a ), (f>± — (h\ 1 hi2), and S^ a (S^ a ) as the pseudospin raising (lowering) operators. It has been 
shown 15 that the constrained electron operator can be mapped exactly using the fermion-spin transformation defined 
with an additional projection operator. However, this projection operator is cumbersome to handle in the actual 
calculations, and we have dropped it in Eq. (2) and in the subsequent calculation. It has also been shown 15 that such 
treatment leads to errors of the order p in counting the number of spin states, which is negligible for small dopings. In 
the two-leg ladder systems, because of the two coupled chains, the energy spectrum has two branches, therefore the 
one-particle spinon and holon Green's functions are matrices, and can be expressed as D(i — j,r — r') = D^{i — j,r — 
t') + a x Dx(i — j,r — t') and g(i — j,T—r') = — j,T — r') + cr x gT(i — j,T — r'), respectively, where the longitudinal 
and transverse parts are defined as D L (i - j,r - r') = -(T T S+(T)Sj~ a (T')), g L (i - j,r - r') = -{T T h ia {T)h) ja {T')) 
and D T (i-j,r- t') = -{T T S+ a {T)Sj a ,{T')), g T (i-j, T - r') = -(T r /i ia (^W')> with a £ a'. Within the charge- 
spin separation, the spin fluctuation couples only to spinons 17 , however, the strong correlation between holons and 
spinons still is included self-consistently through the holon's order parameters entering the spinon's propagator, then 
both holons and spinons are responsible for the spin dynamics. In this case, the spin dynamics of the doped square 
lattice antiferromagnet have been discussed 18 , and the results are consistent with the experiments 9,10 . Following their 
discussions 18 , the dynamical spin structure factor of the doped two-leg ladder antiferromagnet is obtained explicitly 
as, 



S(k, lu) = -2[1 + n B (w)] [2ImL> L (k, u) + 2ImL> T (k, lo)] 

4[l + n B H](4 1) ) 2 ImE^(k,c) 



(4 1} ) 2 - B^ReSW (k, cvW + [BflrnSW (k )W )] 5 



(3) 



where the full spinon Green's function, D^fcu) = ^"'(k.w) - E< s )(k,w), with the mean-field spinon Green's 
function, 
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v = 1,2, ImE^(k,w) = ImE^k,^) + ImS^ s) (k,w), ReE^(k,w) = ReS£°(k,w) + ReE^ s) (k,w), while 
ImE^(k, uj) (ImE^(k, uj)) and ReS^(k, uj) (ReE^(k, uj)) are the imaginary and real parts of the second or- 
der longitudinal (transverse) spinon self-energy, respectively obtained from the holon bubble as (k, uj) = 

(l/i) 2 E ¥ E w v» n -v«(k,P,P» and E^(k,c) = (1/L) 2 E PP < £,w (-l) v+v ' +v " +1 ^u" (k, p, p', w) with L 
is the number of rungs, and 

IW^(k,p,p» - (2t||[cos(p'+p + fc) +cos(p' -fc)] +t ± [(-l)' /+y " + (-1)"+""; 

<2 f f£Uk,P,p') F£)„„(k,p,p') \, (5) 
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where B^ = B k - J ±cB [ X ± + 2xl(-l) v ][e ± + (-1)1, B k = A[(2e|| X jj + X||)cosfc - (e||X|| + 2 X f)], A - 4J|| eff , e„ = 
1 + 2t||0||/J|| eff , ej_ = 1 + 4i_L^j_/ Jj_ e ff , and 

F$ v „(k,p,p') = n F (^)[l - n F (£^)] - n B ) [n F (#> ) - n F (^ + 'J,)], 

*£Uk,p,p') = M^M 1 n F (^)\ + [1 + n B (J k ^)][n F (^) - n F (^)], (6) 

n F (^[ 1/ ' 1 ) and ns(w^) are the fermion and boson distribution functions, respectively the mean-field holon excitations 
= 4t||X||Cosfc + ^ + 2x±ij_(-l)' y+1 , and the mean-field spinon excitations, (w^) 2 = uj\ + A 2 (-l) ,y+1 with 
cj 2 = ^4i(cos/c) 2 + A 2 cosk + A 3 , A 2 = Xicosfc + X 2 , where 

A, = a e||A 2 (x||/2 + £||x|), 

A 2 = -6||A 2 Ke|| X ||/2 + xf)/2 - a(q + Cy/2) - (1 - a)/4] 

- aAJ ±off [e||(Ci + xl) + e±(Cj_ + e X ±)/2], 
A 3 = \ 2 [a(q + efq/2) + (1 - a)(l + ef)/8 - ae||( X ||/2 + e|| X f)/2] 

+ aAJ ±cff [e||e ± Cj_ + 2CJ] + Jj_ eff (ei + l)/4, 

= aAJLcff [(e ± X|| + e||X±)/2 + e||£_L(xl + xf )], 
X 2 = - a AJ ±cff [e||e ± X||/2 + e±(xf + C x ± ) + e\\C±/2] - e ± Jl cS /2, (7) 



with the spinon correlation functions X || = (Sai S ai+fj)> X\\ = ( S ai S i+fj), X± = (SuS 2i ), xl = (SuSZi), C \\ = 
(mEri(S+ +fj S- + -,), q - (l/4)E I>? (^ + ,^ + ,,), C ± = (1/2)E,(^ 2 + ^h + ,), and ^ = (1/2) E^Sf^}. 

In order to satisfy the sum rule for the correlation function (S^S^A = 1/2 in the absence of the antiferromagnetic 
long-range-order, a decoupling parameter a has been introduced in the mean-field calculation, which can be regarded 
as the vertex correction 19 . All the mean-field order parameters, decoupling parameter, and chemical potential have 
been determined self-consistently, as done in the two-dimensional case 18 ' 19 . 

In Fig. 1, we represent the dynamical spin structure factor 5(k, ui) in the (k,w) plane at doping p = 0.16 with 
temperature T = 0.05J|| for t\\/ J\\ — 2.5, (a) t±/t\\ — 0.5, J±/J\\ — 0.5, and (b) t±/t\\ = 0.7, J±/J\\ = 0.7, hereafter we 
use the units of [2ir] . Obviously, an interchain coupling dependence of the incommensurate-commensurate transition 
occurs. To check this point explicitly, the calculated dynamical spin structure factor spectrum has been used to 
extract the doping and interchain coupling dependence of the incommensurability 8, defined as the deviation of the 
peak position from the antiferromagnetic wave vector Q = [1/2], and the results in T = 0.05J|| and to = 0.4Jm for 
t\\/ J\\ = 2.5, t_\_/t\\ = and J±/J\\ = (dash-dotted line), t±/t\\ = 0.3 and J±/J\\ = 0.3 (solid line), t±/t\\ = 0.5 
and J±/J\\ = 0.5 (long-dashed line), fj_/f|| = 0.65 and J±/J\\ = 0.65 (short-dashed line), and t±/tu = 0.69 and 
J±/J\\ = 0.69 (dotted line) are shown in Fig. 2. We therefore find that there are a critical values of t±/t\\ = 0.7 and 
J±/J\\ = 0.7 in T = 0.05J|| and u) = 0.4J||. In the weak coupling regime with t±/t\\ <^ 0.7 and J±/J\\ <C 0.7, the 
commensurate scattering peak near the half-filling is split into two peaks at [(1 ± 8)/2], where 8 increases initially 
with the hole concentration in lower dopings, but it saturates at higher dopings. In this case, spinons and holons 
are more likely to move along the legs of the ladders, rendering the materials quasi-one-dimension. However, the 
range of the incommensurate spin correlation decreases with increasing the strength of the interchain coupling, and 
the commensurate spin fluctuation appears in the whole doping range in the strong coupling regime with tj_/tu > 
0.7 and J±/J\\ > 0.7. Many experimental analyses 1 ' 5 ' 7 have indicated that J±/J\\ > 0.7 for the doped two-leg 
ladder antifcrromagnet Sri4_ x Ca x Cu2404i, this is why that the incommensurate spin correlation is not observed 
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in Sri4_ a; Ca a ;Cu2404i from the inelastic neutron scattering experiments 1,5 ' 7 . For a better understanding of the 
evolution of the commensurate scattering peak with energy in Sri4_ x Ca x Cu2404i, we have made a series of scans 
for the dynamical spin structure factor S*(k, to) at different dopings, and the result for p = 0.2 with T — 0.05 J | for 
t\\/J\\ = 2.5, t±/tii = 0.77 and J±/J\\ = 0.77 is plotted in Fig. 3 in comparison with the experimental data 7 taken 
on Sri4_ x Ca a; Cu2404i (inset) with x = 11.5 (p ~ 0.2) and J±/J\\ = 0.77. This result shows that the commensurate 
spin fluctuation is energy dependent, and the commensurate scattering peak, which is similar to the resonance peak 
in the bilayer cuprate in the normal state 10 , is located at energy lo = 0.93J||. This reflects that the anticipated 
spin gap A s = 0.93J|| w 83.7mcV (Rcf. 7 ) is larger than the spin gap « 32.1meV observed 7 in Sri4_ x Ca x Cu2404i, 
which may mean that the simplest two-leg t-J ladder can not be regarded as the complete model for the quantitative 
comparison with the doped two- leg ladder antiferromagnet. Furthermore, we have also made a series of scans for the 
dynamical spin structure factor at different temperatures, and found that the weight of the peak is suppressed severely 




FIG. 1. The dynamical spin structure factor in the (k,co) plane at p — 0.16 with T — 0.05J|| for t\\/J\\ = 2.5, (a) t±/t\\ = 0.5 
and J±/J\\ = 0.5, and (b) t±/t\\ = 0.70 and J±/J\\ = 0.70. 
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FIG. 2. The doping dependence of incommensurability 8(x) in T = 0.05J|| and lo — 0.4J|| for t\\/J\\ = 2.5, t±/t\\ = and 
J±/J\\ = (dash-dotted line), t±/t\\ = 0.3 and J±/J\\ = 0.3 (solid line), t±/t\\ = 0.5 and J±/J\\ = 0.5 (long-dashed line), 
£_l /i | = 0.65 and J±/J\\ = 0.65 (short-dashed line), and t±/t\\ = 0.69 and J±/J\\ = 0.69 (dotted line). 
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with increasing temperatures. Our present results are in qualitative agreement with the experimental data 7 . On the 
other hand, our present conclusion in the strong coupling regime contradict the numerical result of the isotropic 
two- leg Hubbard ladder 20 , where the equal-time spin structure factor has been calculated using the density-matrix 
renormalization group method, and the result shows that the commensurate peak at the half-filling is split into two 
peaks with dopings. But we want to stress that these results are obtained within the two-leg Hubbard ladder, where 
the Hubbard U is finite. However, in the t-J model, the doubly occupied Hilbert space has been pushed to infinity 
as the Hubbard U — > oo, and therefore the dynamical spin structure factor in the t-J model only describes the lower 
Hubbard band. This may lead to some different results between the Hubbard and t-J models. Of cause, this has to 
be checked by the further density-matrix renormalization group study for the t-J ladder. 

One of the most important features of the spin dynamics in the doped two-leg ladder antiferromagnet is the spin- 
lattice relaxation time 1/Ti 8 . This spin-lattice relaxation time is closely related with the dynamical spin structure 
factor, and can be expressed as, 

1/Tx = 2K B T/(g^%h) lim l(1/JV) £ F Q 2 (k) X "(k, (8) 

k 

with g is the g factor, jiQ is the Bohr magneton, F a (k) is the form factors, and dynamical spin susceptibility x"(k, oj) is 
related to the dynamical spin structure factor by the fluctuation-dissipation theorem as, x"(k, w) = (1 — e~ t3uj )S(k, u>). 
The form factors have dimension of energy, and magnitude determined by atomic physics, and k dependence deter- 
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FIG. 3. The dynamical spin structure factor at the AF wave vector Q = [1/2] for p — 0.2 with T = 0.05J|j for tu/Ju = 2.5, 
t±/t\\ — 0.77 and J±/J\\ = 0.77. Inset: the experimental result 7 on Sri4_ I Ca I Cu2404i with x = 11.5 (p ft* 0.2) and 
J±/J\\ = 0.77 taken from Ref. 7 . 
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FIG. 4. The temperature dependence of spin-lattice relaxation time 1/Ti in both logarithmic scales at p = 0.20 for 
t\\/J\\ — 2.5, t±/t\\ — 0.77 and J±/J\\ = 0.77. Inset: the experimental result on Sri4- :E Ca :E Cu2404i with x = 11.5 (p ~ 0.2) and 
J±/J\\ = 0.77 taken from Ref. 8 . 
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mined by geometry. Since the strong short-range commensurate spin fluctuation in the strong coupling regime in the 
dynamical spin structure factor (3), the main contribution to 1/T\ comes from the region around the antiferromag- 
netic wave vector, therefore we can set F a (h) as constant without loss the generality 13 . In this case, the spin-lattice 
relaxation time 1/T\ has been evaluated and the result at p — 0.20 for t\\/ J\\ = 2.5, t±/t^ = 0.77 and J±/J\\ = 0.77 
is plotted in Fig. 4 in comparison with the experimental data 8 taken on Sri4_ :c Ca a ;Cu2404i (inset) with x = 11.5 
(p « 0.2) and J±/J\\ = 0.77, where we have chosen units h = Kb = 1. This result indicates that in low temper- 
atures the spin-lattice relaxation time decreases exponentially with decreasing temperatures, in agreement with the 
experiments 8 . 

The dynamical spin structure factor in Eq. (3) has a well-defined resonance character, where S(k, lo) exhibits peaks 
when the incoming neutron energy lo is equal to the renormalized spin excitation E\ = (lo^) 2 + B^ReS^,(k, E^), 
i.e., W(k Cl uj) = [lo 2 - (lo^) 2 - B^Re£^(k c , lo)} 2 = (lo 2 - E 2 ^) 2 ~ for certain critical wave vectors k c . Then 

the height of these peaks is determined by the imaginary part of the spinon self-energy 1/Im£^(k c ,w). This 
renormalized spin excitation is doping, energy, and interchain coupling dependent. In the present spinon self- 
energy Re£^(k,w) = ReS^ } (k,w) +ReS^ s) (k,w), ReE^ s) (k,w) < favors the one-dimensional behaviors, while 
RcSy ' (k, lu) > characterizes the quantum interference between the chains in the ladders, therefore there is a com- 
petition between ReE^ (k, lo) and ReS^ (k, lo) . In the weak coupling regime, the main contribution for ReS^(k, lo) 

may come from ReS^ (k, lo), and spinons and holons arc more likely to move along the legs, then the incommensurate 
spin correlation emerges, where the essential physics is almost the same as in the two-dimensional t-J model 18 . Near 
the half-filling, the spin excitations are centered around the antifcrromagnetic wave vector [1/2], so the commensurate 
antifcrromagnetic peak appears there. Upon doping, the holes disturb the antiferromagnetic background. Within the 
fermion-spin framework, as a result of self-consistent motion of holons and spinons, the incommensurate spin corre- 
lation is developed beyond certain critical doping, which means, the low-energy spin excitations drift away from the 
antiferromagnetic wave vector, where the physics is dominated by the spinon self-energy ReE^ (k, lo) renormalization 
due to holons. However, the quantum interference effect between the chains manifests itself by the interchain coupling, 
i.e., this quantum interference increases with increasing the strength of the interchain coupling. Thus in the strong 
coupling regime, ReS^ (k, lo) may cancel the most incommensurate spin correlation contributions from ReS^ (k, lo), 
then the commensurate spin fluctuation appears. In this sense, the interchain coupling is a crucial role to determine 
the symmetry of the spin fluctuation in the doped two- leg ladder antiferromagnet. Since the height of the peaks is 
determined by damping, it is fully understandable that the weight of the peak is suppressed as the temperature are 
increased. 

To conclude we have discussed the spin dynamics of the doped two- leg ladder antiferromagnet within the t-J 
ladder. Our calculations clearly show a crossover from the incommensurate spin correlation to commensurate spin 
fluctuation characterized by the rung to chain coupling. However, in the regions of parameter space given by experi- 
ments, the t-J ladder can correctly reproduce all main magnetic features of the doped two-leg ladder antiferromagnet 
Sri4_ :c Ca^Cu2404i, including the energy dependence of the neutron scattering peak position and exponential decrease 
of the nuclear spin-lattice relaxation rate. 
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